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Abstract 

A basis set expansion is performed to find the eigenvalues and wave functions for an electron on 
a toroidal surface T 2 subject to a constant magnetic field in an arbitrary direction. The evolution 
of several low-lying states as a function of field strength and field orientation is reported, and a 
procedure to extend the results to include two-body Coulomb matrix elements on T 2 is presented. 
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1. INTRODUCTION 



Quantum dots with novel geometries have spurred considerable experimental and the- 
oretical interest because of their potential applications to nanoscience. Ring and toroidal 
structures in particular have been the focus of substantial effort because their topology 
makes it possible to explore Ahranov-Bohm and interesting transport phenomena Q ill J. 
Toroidal InGaAs devices have been fabricated [ESQ Bl and modelled, and toroidal car- 
bon nanotube structures studied by several groups 

This work is concerned with the evolution of one-electron wave functions on T 2 in response 
to a static magnetic field in an arbitrary direction. The problem of toroidal states in a 
mag netic field has been studied with various levels of mathematical sophistication. Onofri 



12j has employed the holomorphic gauge to study Landau levels on a torus defined by a strip 
with appropriate boundary conditions and Narnhofer has analyzed the same in the context 
of Weyl algebras Here, the aim is to do the problem with standard methodology: 

develop a Schrodinger equation inclusive of surface curvature, evaluate the vector potential 
on that surface, and proceed to diagonalize the resulting Hamiltonian matrix. 

As noted in [3], ideally one would like to solve the A-electron case, but the single 
particle problem is generally an important first step, and while the N electron system on 



flat and spherical surfaces has been studied [1^ lla llll lla llSl; 20], the torus presents its own 



difficulties. In an effort to partially address this issue, the evaluation of Coulombic matrix 
elements on T 2 is also discussed here. 

This paper is organized as follows: in section 2 the Schrodinger equation for an electron 
on a toroidal surface in the presence of a static magnetic field is derived. In section 3 a 
brief exposition on the basis set employed to generate observables is presented. Section 4 
gives results. Section 5 develops the scheme by which this work can be extended to the two 
electron problem on T 2 , and section 6 is reserved for conclusions. 



2. FORMALISM 

The geometry of a toroidal surface of major radius R and minor radius a may be param- 
eterized by 

r(0,0) = W(6)p + a sm6k (1) 

with 

W = R + a cos6, (2) 
p = cos^i + sin^j. (3) 

The differential of Eq.(l) 

dr = ad6 6 + Wd<j><j> (4) 
with 6 = —sin6p + cos^k yields for the metric elements gij on T 2 

9ee = a 2 (5) 



2 



= W 2 . (6) 
The integration measure and surface gradient that follow from Eqs. (5) and (6) become 

^gdq l dq 2 -> aWdOdcj) (7) 

and 

„1 d A 1 d 

add W o(p 

The Schrodinger equation with the minimal prescription for inclusion of a vector potential 
A is 

1 ( H \ 2 
H = — ( -V + gA | * = EV. (9) 

2m \i J 

The magnetic field under consideration will take the form 

B = Bii + B k, (10) 

which by symmetry comprises the general case. In the Coulomb gauge the vector potential 
A(9, 0) = ~B x r expressed in surface variables reduces to 

A(6,<f>) = -[Si(Wsin0cos0 + a sin 2 flsin</>)0 + {B W - B x a sin^cos0)]0 



+Bi (Fsin0 sin 9 — a cos^sin^ sin0)n. 



(11) 



with n = cf) x 0. The normal component of A contributes a quadratic term to the Hamilto- 
nian but leads to no differentiations in the coordinate normal to the surface as per Eq.(8). 

e is constrained to 



There is a wealth of literature concerning curvature effects when apartic 
a two-dimensional surface in three-space 21, 22,|i3, 24, 25, 26, 27, 28 L 



34, 



35, 



36, 



37j,|38j, including some dealing with the torus specifically |39j, but the scope of 
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33, 



this work will remain restricted to study of the Hamiltonian given by Eq. (9). 

The Schrodinger equation (spin splitting will be neglected throughout this work) is more 
simply expressed by first defining 

a = a/R 
F = 1 + a cos9 



7o 
7i 



B ttR 2 
B^R 2 



In = — 
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2mEa 2 



after which Eq. (9) may be written 

' d 2 a sin 9 d a 2 d 2 
£P9 F~d9 + F2 



+ z I r a -p- sm0cos0 I 
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.9 

+zarisin0(a + cos6>) — 
<76> 



r 2 a 2 F 2 T(a 2 F 2 , ± a 1 sin 2 0\ r r ia a F . „ , 
1 sin H — — 1 H sm0cos0 



F 2 



(12) 
(13) 



3. CALCULATIONAL SCHEME 

To proceed with a basis set expansion, Gram-Schmidt (GS) functions orthogonal over 
the integration measure F = 1 + a cos9 must be generated. Fortunately, it is possible 
to construct such functions almost trivially. The method for doing so has been described 
elsewhere j4o| , so only the salient results will be presented below. 

The 7~i = 0, 9 — > —9 invariance of H T suggests that the solutions of the Schrodinger 
equation be split into even and odd functions, and the primitive basis set can be taken to 
possess this property; 

u n {9) = -^=cos[n9], v n (9) = -^-=sm[n9]. (14) 

'7T \/7l 



The GS functions will take the form 



v m {9) 



with the cxm given by (momentarily supressing the parity superscripts [4l| ) 

c Km = (-) K+m N K (N K _ 1 (3N K _ 1 )(N K . 2 (3N K _ 2 )...(N m (3N m ) (16) 
and the normalization factors Nk determined from 

= T^m (17) 

starting from N = y 1/2 for positive parity states and N\ — 1 for negative parity states. 
The Kv th basis state is attained by appending azimuthal eigenfunctions onto the GS func- 
tions described above, 
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The matrix 



tKKvv 



> 



(19) 



is then easily constructed since the matrix elements can all be written in closed form (see the 
Appendix), and the eigenvalues and eigenvectors determined here with a 30 state expansion 
for each ^-parity. The ordering convention adopted for the states was taken as 



4. RESULTS 

Rather than present a large number of tables conveying little useful information per unit 
page length, the focus will be on indicating how some low-lying states evolve as a function 
of magnetic field strength for two distinct orientations. Some remarks will also be made 
regarding the general trend seen for higher excited states. Here the ratio a = a/R was set 
to 1/2 as a compromise between smaller a where the states tend towards decoupled ring 
functions and larger a which are less likely to be physically realistic. 

Fig. 1 illustrates the evolution of the energy eigenvalue for five low-lying states as a 
function of r with T\ = 0. The states are all distinct and are labelled in the caption. Not 
shown are values trivially obtained from the ±z/-Bo splitting arising from Bq = 0, v ^ 
degeneracy. It is interesting that level crossings with attendant movement towards a ground 
state with different Kv occurs near integer values of To, though it is not immediately clear if 
this is of real significance. It is also of interest to show the sensitivity of the dependence of 
\J/*\J/.F on field strength. Fig. 2 shows that even for moderate field values (r = 5 corresponds 
to a field of 1.3 T for a torus with R = 50 nm) the large effective flux as compared to atomic 
or molecular dimensions causes substantial modification to \I/*\I/F in the ground state. 

The results given in Figs. 1-2 were for a field configuration that did not mix azimuthal 
basis states. To investigate an asymmetric case, let r = and vary T\ wherein no field 
threads the torus. Fig. 3 is analogous to Fig. 1 as described above with the notable 
exception that the v splitting is non-trivial (hence fewer distinct states are shown), and 
level crossings occur at larger values. There is no plot analogous to Fig. 2 for larger values 
of Ti with t = because dependence increases rapidly in the eigenstates even for small 
Ti] a contour plot is preferable. Figs. 4 and 5 show contour plot results for two states at 
three field strengths. Note that there is slightly more dependence in 9 when T\ = for the 
state displayed in Fig. 6 than in Fig. 5; the integration measure acts to cancel the angular 
variation of the state displayed in Fig. 5. 
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yielding a Hamiltonian matrix blocked schematically into 
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5. EXTENSION TO COULOMB INTEGRALS 

The two-electron problem on T 2 is complicated by the inability (at least by the author) 
to find a transformation that decouples the relative electron motion from their center of 
mass motion as is easily done on R 2 0|. The obvious transformations do not lead to any 
advantage over the method adopted by workers long used in atomic and molecular physics, 
which is to evaluate the two-body matrix elements (supressing spin indices and physical 
constants) 

J J $*(r 1 )$*(r 2 )\/(r 1 ,r 2 )(l-P 12 )$ fe (r 1 )<l> i (r 2 )rf 3 r 1 rf 3 r 2 (20) 

with 

V{v u v 2 ) = 4tt]T — ^^YLMWiMYZufoM- (21) 

L,M T> 

Eq. (20) can be adopted on T 2 subject to some peculiarities which are due to the restriction 
of r 1; r 2 to a surface. Eq. (20) on T 2 with the notation employed in section 3 becomes 

... r VpuyWi, 4>i)%v 2 (02, <hMr 1} r 2 )(l - P 12 ) 
Jo Jo 

^ Rvz {e 1 ,4>i)^s^e2A2)F{e 1 )F{e 2 )de 1 de 2 d ( j> 1 d ( j> 2 . (22) 

Consider the direct term; in terms of a spherical coordinate system centered at the middle 
of the torus (r sin# s cos0, r sin# s sin0, r cos# s ) 

r sin^g = R + a cos9 (23) 



r cos^s = a sin^ (24) 



r = VR 2 + a 2 + 2aR cos# (25) 

and defining 

pij^PMMo) (26) 

gives for Eq. (22) after some manipulation 

p LM (e sl (e 1 ))p LM (e S2 mmojfm { ^f^^J^ de ^- ( 27 ) 

The arguments of the Plm are evaluated with 

/ R + a cos9i \ . . 

"•■= arcta,i K^- • (28) 
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To evaluate the integral care must be taken with the >, < character of the radial factor in 
the integrand. One way to proceed is as follows: 

1. Fix 61 = 0. At this point r\ is at its maximum; integrate the integrand of Eq.(27) 
numerically over d9 2 from [0, 27r] by some suitable method to attain a value labelled by, say, 
Go(0i). 

2. Set 9\ = 5. Integrate d02 with r> = r±, r< = r 2 over the interval [5, 2n — 8], then set 
r> = r 2 and r< = r\ from [2ti — 5,5]. This is G$(9\). 

3. Repeat the second step until the entire interval around the toroidal cross section is 
covered. A table [G (9x), Gs{9\), G 2 s(9i)...] results that can then be integrated numerically. 
The exchange term proceeds similarly; only the densities need modification. 



6. CONCLUSIONS 

This work presents a method to calculate the spectrum and wave functions for an electron 
on T 2 in an arbitrary static magnetic field. Aside from the character of the solutions and 
numerical data, perhaps the main result of this paper has to do with the ease with which 
an arbitrarily large number of GS states can be trivially generated. Because every physical 
interaction can eventually be expressed as a periodic function on T 2 , matrix elements for 
the interaction may then be evaluatued in closed form; hence the only restriction to doing 
any problem on T 2 is matrix inversion. 

The procedure employed here to generate observables lends itself to easy incorporation of 
an arbitrary number of surface delta functions or other type of potential. This is important 
because on a closed nanotube, in contrast to a macroscopic crystal, there are a relatively 
small number charge carriers, so the continuum approximation may break down for smaller 
torii. Clearly, the magnetic field treated here will not be sufficient to comprise the general 
case as soon as any sort of lattice structure breaking azimuthal symmetry is imposed on the 
torus. However, the extension is simple to implement in Eq. (12), requiring only a few more 
terms. It would be also be interesting to see the extension of the static case discussed here 
to a time dependent laser control problem of the type in 0. 

Some remarks should be made regarding the curvature potential Vc well known to workers 
in the field of quantum mechanics on curved surfaces. It was shown in (44| that a full three 
dimensional treatment of the problem of a particle near, but not necessarily restricted to 
T 2 , yields a spectra consistent with inclusion of Vc added to the two dimensional surface 
Hamiltonian. Here the potential could not be included without substantially increasing the 
scope and complexity of the problem undertaken. It was shown in [45] that the inclusion 
of a vector potential precludes a separation of variables into surface and normal degrees of 
freedom; A added to the Schrodinger equation requires solving coupled differential equations 
in the surface and normal variables, or if a basis set expansion is employed, a much more 
complicated procedure to generate three-dimensional GS states. 
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Appendix 

This appendix gives closed form expressions for the matrix elements needed to construct 
the matrix = LK^v\H t \K^v). First let 

tKKvv \ I ' I / 

3 , 

Pi = 1 + -a 2 

3 o 

P 2 = 3a + +-a s 
4 



/(«) = 



v 7 !^ 2 - l 



and define 

Each operator in Eq. (9) will connect either only like parity states or opposite parity states; 
no single operator will do both. The matrix elements that connect like positive parities are 



d 2 



d 2 9 



K K 



K + u) = 7i ^ ^Cx m c^ n (-n 2 )Ap_^[(A m+n - A m _ n ) + 

/ m =0 n=0 

~(A m+n+ i + A m _ n+1 + A m+n _! + A m _ n _i)] 



(Al) 



a . d 
F Sm9 W 



CKm c Kn Tl Ap_ I/ ( A m+n _i + A m _ n+ i — A m _ n _i) 

(A2) 



m=0 n=0 



a 2 d 2 



F 2 d 2 <p 



K + A = -=== J2J2 c ^cku (-^ 2 ) A^ u [r +m (a) + p n - m Ka)} (A3) 
I vi- at m=0n=0 



K + i> 



d 



tT a 



K K 



K + v ) = 7ir a 2 ^ ^2 c R m CKn{-v)A p _ u [{A m+n + A m _ n )+ 



m=0 n=0 



a 



(A m+n+ i + A m _ n+ i + A m+n _i + A m _ n _i)] 



(A4) 



T 2 f> 2 



•y •) K K 



K + V ) = °— J2 Yl C R m CKnA V - v 



m=0 n=0 
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[P 1 A m _ n + -(P 2 (A m+n _ 1 + A m _ n+1 + A m _ n _x)+ 
^3(A m+n _2 + A m _ n+2 + A m _„_ 2 ) 

^m— n+3 

+ A m _ n _ 3 ))] 



(A5) 



Ti 2 q; 2 ~ 9 



1 

2' 



9 _ftT ii - 

cr v ^ x - 

2 , 2 / C Km c Kn 




4- 


m=0 n=0 






^m+n— 1 "I - A m _ n _|_x 


A m _ n _ 


-0+ 


> + A m _ n+2 + A m _ 


n-2) 




+ A m _ n+ 3 + A m _ 


-3))] 





(A6) 



r l 2 « 4 • 2/1 

— sin^ 



nT 2 a 4 R K 

K~v) = 1 — ^2^2c Rm c Kn A p _ u [A m+n + A m _ n + 
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m=0 ra=0 



(A m+n+ i + A m _ n+ i + A n _ m+ i + Ai_ 



+ A m _ n+2 + A n _ m+2 + A 2 - m _„) 



a 



-j-(A m+n _ 3 + A m _ n+3 + A m _ n _ 3 + A 3 _ m _ n ))] . 



(A7) 



The negative to negative terms are 



d 2 



d 2 e 



K K 



K v \ = 7T ^ ^2 d KmdKn(-n 2 )A^ u [A m _ n + ^(A m _ n+ i + A m _ n _i)] (A8) 



m=l n=l 



a . d 
F Sm9 d6 



v K K 

K~v \ = a- ^2^2d Rm d Kn n A p _„(A m _ n+1 - A n _ m+1 ) (A9) 

I ■KK) 1 t~l 1 



m=l n=l 



a 2 d 2 



F 2 d 2 (j) 



9 K K 

a 71 



VT^ 2 



J2Y, d *rndKn (^ 2 ) A P ^[P +m (a) - /l»-™l(a)] (A10) 



m=l n=l 



if"*/ 
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2 



K u) = 7TT a 2 ^ ^ dK m d K n{-v)&v-v [A m - n + 

' m=l n=l 

(A m+n+ i + A m _ n+ i + A m+n _i + A m _ n _i)] 



(All) 
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1 



9 9 K K 



K v \ = ^— ^ d Rm d Kn ^u-u 



m=l n=l 



n(A(A m -H-l + A n _ m _! — A!_ m _ n ) + 

P 3 (A m— n— 1 A n _ m _2 A2— m— rt) 
+-P4(A m _ n _i + A n _TO-2 _ A 3 _ n _ m ))] 



(A12) 



T l « 2 r.2 • 2 



9 o K K 
7TT{a 



1 

2' 



m=l 


n=l 




m—n—1 


A ra _ m _i 


Ai_ m _ n 


+ A n _ 


-m-2 — A 2 - 


m— rt) 




-m-3 + A 3 _ 


m— n))] 



{Aid) 



-2„A 



T l a -2/1 

— sin^ 



7rr 2 Q ,4 ^ ^ 

/V "? ' ) = | — ^2 ^ d Rm d Kn A p _ u [A m _ n + 

m=l n=l 



■^(A m _ n+ i + A n _ m+1 — A!_ m _ n ) + -(— A m _ n+2 — A m _ n _ 2 + A 2 _- m _ n ) 



Oi 

+-r(— A m _ n+3 — A m _ n _ 3 + A 3 _ m _ n ))]. 



The matrix elements connecting different parities are 



(AU) 



K + 



v 



— i— — sm.9cos6— - 
F d<j) 



K \ = nnQ_ ^2 ^ C Rm d K n{-v){&v-v+l + ^u-u-l) 



K K 



m=0 n=l 



[A m+n+ i + A n _ m+i — A m _ n+1 — A m+nl _!] 



(A15) 



K + 



v 



d 

iarisin0(o; + cos#) — 
06 



K K 



K v) = t x o?tx ^ ^ c Rm d K n{n)(/± v - V+1 - A„_p_i) 



m=0 ra=l 



r 3a; (1 + a 2 ) 



+ |(A m +n+2 + A m _ n _ 2 + A m+n _ 2 + A m _„_ 2 )] 



(A16) 



ff 



K + i 



T Ticr 



sin^cos0F 



K~v ) = 



TqT\OL 7T 



m=0 n=l 



la 2 

+ T)(A m + n -l — A m+n+1 + A m _ n _! — A m _ n+1 ) + 



— (A m+n _2 — A m+n+2 + A m _„_ 2 — A m _ n+2 ) 



+-^-(A n+m+ i — A m+n _i + A m _ Tt+ i — A m _ n _i 



+A^ 



+rra— 3 A m _|_ n _|_3 -\~ A n _ m _3 



A n _ m+ 3)] . 



0417) 



The negative to positive elements are obtained by interchanging all indices, or equivalently, 
by taking their transpose. 
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Figure Captions 



Fig. 1: e as a function of To for five low-lying states. Diamonds correspond to the \vK > 
= |00 + > state, stars to | — 10 + >, squares to | — 20 + >, triangles to 1 01 > and circles to 
|01+ >. 

Fig. 2: Evolution of ^/*^/F for the |00 + > state given for r = (thin line), r = 2.5 
(medium line) and t = 5.0 (thickest line). 

Fig. 3: e as a function of T\ for five low-lying states. Diamonds correspond to the \uK ± > 
= |00 + > state, stars/squares to | — 10 + >, and triangles/circles to | — 20 + >. 

Fig. 4: Sequential evolution of the ^*^F = F\00 + > state on T 2 for r = 0, r = 2.5 and 
r = 5.0. 

Fig. 5: Sequential evolution of the ty*tyF — F\ - 10 + > state on T 2 for r x = 0, t x = 2.5 and 
Ti = 5.0. The ground state variation in 9 is partially cancelled by the integration measure. 
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Fig. 3 
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